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@ Computertomographie-Verfahren mit kegelformigem Strahlenbiindel 

(57) Die Erfindung betrifft ein Computertomographie-Ver- 
fahren, bei dem der Untersuchungsbereich auf einer 
Kreisbahn mit einem kegelformigen Strahlenbundel ab- 
getastet wird. Der Rechenaufwand fu r die Rekonstruktion 
der Absorptionsverteilung im Untersuchungsbereich laftt 
sich dabei durch folgende Schritte vereinfachen: 

a) Rebinning der Mefcdaten zu einer Anzahl von Gruppen, 
wobei jede Gruppe mehrere zueinander und zur Rotati- 
onsachse parallele Ebenen umfaftt, in denen sich je ein 
Strahlenfacher befindet, 

b) eindimensionale Filterung der durch das Rebinning er- 
zeugten Daten einer jeden Gruppe in Richtung senkrecht 
zur Richtung der Ebenen, 

c) Rekonstruktion der raumlichen Verteilung der Absorp- 
tion aus den gefilterten Daten verschiedener Gruppen. 
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Beschreibung 

Die Erfindung betrifft ein Computertomographie-Verfah- 
ren mit den Schritten 

5 

- Erzeugen eines kegelformigen, einen Untersu- 
chungsbcrcich bzw. cin darin bcfindlichcs Objckt 
durchsetzenden Strahlenbundels mit einer Strahlen- 
quelle, 

- Erzeugung einer kreisformigen, eine Rotation um 10 
eine Rotations achse umfassenden Relativbewegung 
zwischen der Strahlenquelle einerseits und dem Unter- 
suchungsbereich bzw. dem Objekt andrerseits, 

- Akquisition von MeBdaten, die von der Intensitat in 
dem Strahlenbiindel jenseits des Untersuchungsberei- 15 
ches abhangen, mit einer Detektoreinheit wahrend der 
Relativbewegung, 

- Rekonstruktion der raumlichen Verteilung der Ab- 
sorption innerhalb des Untersuchungsbereiches aus 
den akquirierten MeBdaten. 20 

AuBerdern bezieht sich die Erfindung auf einen Compu- 
tertomographen zur Durchfiihrung dieses Verfahrens. 

Ein solches Verfahren ist aus dem Aufsatz von L. A. Feld- 
kamp et al "Practical cone-beam algorithm", Journal of Op- 25 
tical Soc. Am. A/VoL 1, No 6, 1984, pp. 612-619 bekannt. 
Das bekannte Verfahren besteht im Prinzip aus den folgen- 
den Schritten: 

a) Fur jede Strahlenquellenpositionen werden alle 30 
MeBwerte mit einem Gewichtungsfaktor multipliziert, 
der dem Kosinus des Winkels entspricht, den der Strahl 
langs dessen der MeBwert erfaBt wurde, mit dem Zen- 
tralstrahl einschlieBt. 

b) Die so gewichteten MeBwerte werden einer Hoch- 35 
paBfilterung unterzogen. 

c) Die MeBwerte werden entlang der Strahlen, langs 
derer sie gemessen worden sind in den Untersuchungs- 
bereich ruckprojiziert. Der Beitrag eines MeBwerts zu 
dem Absorptionswert eines Voxels muB dabei mit ei- 40 
nem Faktor gewichtet werden, der vom Abstand des 
betreffenden Voxels von der Strahlenquellenposition 
abhangt. 

Da dieser letzte Schritt fur samtliche Voxel des zu rekon- 45 
struierenden Volumens und fiir samtliche Strahlenquellen- 
positionen durchgefiihrt werden muB, benotigt er viel Re- 
chenzeit. Ein weiterer Nachteil des bekannten Verfahrens 
ist, daB man damit die Absorption nur in denjenigen Voxeln 
rekonstruieren kann, die wahrend der Untersuchung unun- 50 
terbrochen von Rontgenstrahlung getroffen worden sind. 
Diese Voxel liegen in einem zur Rotations achse konzentri- 
schen diskusformigen Bereich. Erwiinscht ist jedoch die Re- 
konstruktion in einem ebenen scheibenformigen Bereich. 
Beschrankt man aber die Rekonstruktion auf den scheiben- 55 
formigen Bereich in dem kegelformigen Strahlenbundel, der 
aus alien Strahlenquellenpositionen von Rontgenstrahlung 
getroffen wird, dann ergibt sich lediglich ein sehr schmaler 
Rckonstruktionsbcrcich . 

Aufgabe der vorliegenden Erfindung ist es ein Computer- 60 
tomographic- Verfahren der eingangs genannten Art so aus- 
zugestalten, daB der Rechenaufwand verringert wird und die 
Rekonstruktion der Absorptions verteilung in einer dickeren 
ebenen Scheibe des Untersuchungsbereichs moglich wird. 
Diese Aufgabe wird ausgehend von einem Verfahren der 65 
eingangs genannten Art erfindungsgemaB durch die folgen- 
den Schritte gelost: 
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a) Rebinning der MeBdaten zu einer Anzahl von Grup- 
pen, wobei jede Gruppe mehrere zueinander und zur 
Rotationsachse parallele Ebenen umfaBt, in denen sich 
je ein Strahlenfacher befindet, 

b) Eindimensionale Filterung der durch das Rebinning 
erzeugten Daten einer jeden Gruppe in Richtung senk- 
rccht zur Richtung der Ebenen, 

c) Rekonstruktion der raumlichen Verteilung der Ab- 
sorption durch Ruckprojektion der gefilterten Daten 
von mehreren Gruppen. 

Wahrend bei dem bekannten Verfahren die in den einzel- 
nen Strahlenquellenpositionen gewonnenen MeBwerte je- 
weils unmittelbar einer Filterung unterzogen werden, erfolgt 
bei dem erfindungsgemaBen Verfahren zunachst ein Rebin- 
ning, bei dem in unterschiedlichen Strahlenquellenpositio- 
nen erfaBte Strahlenfacher gruppenweise zusammengefaBt 
und erst danach gefiltert werden. Der Vorteil des zusatzli- 
chen Rebinning besteht darin, daB bei der anschlieBenden 
Rekonstruktion die Absorptionsverteilung die Beitrage der 
einzelnen MeBwerte zu einem Voxel nicht mit einem ab- 
standsabhangigen Faktor multipliziert werden miissen. Da- 
durch vereinfacht sich das Rekonstruktionsverfahren ganz 
deutlich. Es ergibt sich auch eine verbesserte Bildqualitat. 
AuBerdern ergibt sich ein giinstiges Verhaltnis zwischen 
dem bestrahlten Volumen und dem Volumen, innerhalb des- 
sen die Rab sorptions verteilung rekonstruiert werden kann. 

Die Erfindung geht weiter von der Erkenntnis aus, daB 
sich die Absorption in jedem Voxel rekonstruieren laBt, 
wenn es iiber einen Winkelbereich von wenigstens 180° be- 
strahlt worden ist. Es laBt sich zeigen, daB alle Voxel, die 
diese Bedingung erfiillen, in einer ebenen, zur Rotations- 
achse senkrechten Scheibe liegen, deren Dicke groBer ist als 
die Abmessungen an den Randern des diskusformigen Be- 
reichs, in dem die Absorption mit dem bekannten Verfahren 
rekonstruiert werden kann. 

Das in Anspruch 2 vorgesehene Rebinning auf einen ebe- 
nen und rechteckigen virtuellen Detektor erleichtert die 
nachfolgende eindimensionale Filterung ganz wesentlich. 

In der Scheibe, in der die Absorptionsverteilung rekon- 
struiert werden kann, befinden sich Voxel, die aus einem 
Winkelbereich von mehr als 180° bestrahlt worden sind. Da 
fiir die Rekonstruktion aber nur ein Bestrahlungs winkelbe- 
reich (das ist der Winkelbereich, den die (Parallel-) Projek- 
tion der Strahlen von der Strahlenquelle zu einem Voxel auf 
eine zur Rotationsachse senkrechte Ebene in dieser Ebene 
bedecken, bzw. den die Komponenten der Vektoren von der 
Strahlenquelle zu dem Voxel in der Rotationsebene der 
Strahlenquelle iiberstreichen) von 180° erforderlich ist, kon- 
nen MeBwerte aus einigen Strahlenquellenpositionen unbe- 
rucksichtigt bleiben. Anspruch 3 beschreibt nun eine einfa- 
che Moglichkeit, wie man die Strahlenquellenpositionen er- 
mitteln kann, deren MeBdaten fiir die Rekonstruktion der 
Absorptionen in einzelnen Voxeln heranzuziehen sind. 

Grundsatzlich kann man die Absorptionsverteilung in der 
erwahnten Scheibe rekonstruieren, indem man fur alle Voxel 
nur so viele Strahlenquellenpositionen zur Rekonstruktion 
heranzieht, daB sich gerade ein Bestrahlungswinkelbereich 
von 180° ergibt. Es gibt in dieser Scheibe jedoch auch Voxel 
mit einem Bestrahlungswinkelbereich von 360°. Wenn man 
diese Voxel ebenfalls nur mit einem Bestrahlungswinkelbe- 
reich von 180° rekonstruieren wiirde, wurde sich fiir diese 
Voxel nicht das bestmogliche Signal/Rauschverhaltnis erge- 
ben. Anspruch 4 beschreibt daher ein hybrides Rekonstruk- 
tionsverfahren, bei dem das Volumen der rekonstruierbaren 
Scheibe in zwei Teilvolumina unterteilt wird: in ein erstes 
Teilvolumen, das nur Voxel umfaBt, die aus alien Strahlen- 
quellenpositionen bestrahlt worden sind und in ein zweites 
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Teilvolumen, bei dem der Bestrahlungswinkelbereich klei- 
ner ist als 360°. Bei der Rekonstruktion der Absorption in 
den Voxeln des ersten Teilvolumens werden samtfiche MeB- 
werte herangezogen und bei der Rekonstruktion der Voxel in 
dem zweiten Teilvolumen nur die MeBwerte aus einem Be- 5 
strahlungswinkelbereich von 180°. 

Anspruch 5 bcschrcibt cincn Computcrtomographcn zur 
Durchfuhrung des erfindungsgemaBen Verfahrens. 

Anspruch 6 beschreibt eine Ausgestaltung des erfin- 
dungsgemaBen Verfahrens, die darauf abstellt, daB ein Teil 10 
der MeBwerte fur die Rekonstruktion nicht herangezogen 
wird. Diese MeBwerte gehoren zu Strahlen, die nicht beno- 
tigt werden, wenn man nur die Ab sorptions verteilung in ei- 
ner ebenen Schicht rekonstruieren will. Diese Strahlen ha- 
ben in einer zur Rotations achse senkrechten Ebene einen 15 
kleinen Offnungswinkel und in einer die Rotations achse 
enthaltenden Ebenen einen relativ groBen Offnungswinkel. 
Wenn man das kegelformige Strahlenbiindel durch die in 
Anspruch 6 angegebene Gestaltung der Kollimatoranord- 
nung so ausformt, daB diese Strahlen gar nicht erst entste- 20 
hen, wird die Strahlenbelastung fiir einen im Untersu- 
chungsbereich befindlichen Patienten verringert. Dieser 
Vorteil ergibt sich auch dann wenn man die Absorptions ver- 
teilung in anderer Weise rekonstruiert, z. B. durch ein ART- 
Verfahren. 25 

Die Erfindung wird nachstehend anhand der Zeichnungen 
naher erlautert. Es zeigen: 

Fig. 1 Einen Computertomographen mit dem die Erfin- 
dung ausfiihrbar ist. 

Fig. 2 Ein Ablaufdiagramm fiir das erfindungsgemaBe 30 
Rekon struktion s verf ahren . 

Fig. 3 Ein in einer Strahlenquellenposition erzeugtes ke- 
gelformiges Strahlenbiindel. 

Fig. 4 Die durch das Rebinning gebildeten Strahlenfacher 
in parallelen Ebenen. 35 

Fig. 5 Einen Querschnitt durch diesen Strahlenfacher. 

Fig. 6 Die geometrischen Verhaltnisse im Untersuchungs- 
bereich. 

Fig. 7 Ein Diagramm zur Auswahl von Strahlenquellen- 
positionen. 40 

Der in Fig. 1 dargestellte Computertomograph umfaBt 
eine Gantry 1, die um eine parallel zur z-Richtung, des in 
Fig. 1 dargestellten Koordinatensystems verlaufende Rotati- 
onsachse 14 rotieren kann. Dazu wird die Gantry von einem 
Motor 2 mit einer vorzugsweise konstanten, aber einstellba- 45 
ren Winkelgeschwindigkeit angetrieben. An der Gantry ist 
eine Strahlenquelle S, beispielsweise ein Rontgenstrahler, 
befestigt. Dieser ist mit einer Kollimatoranordnung 3 verse- 
hen, die aus der von der Strahlenquelle S erzeugten Strah- 
lung ein kegelformiges Strahlenbundel 4 ausblendet, d. h. 50 
ein Strahlenbundel, das sowohl in z-Richtung als auch in ei- 
ner dazu senkrechten Richtung (d. h. in einer zur Rotations- 
achse senkrechten Ebene) eine von Null verschiedene, end- 
liche Ausdehnung hat. 

Das Strahlenbundel 4 durchdringt einen Untersuchungs- 55 
bereich 13, in dem sich ein Patient auf einem Patientenlage- 
rungstisch (beides nicht naher dargestellt) befinden kann. 
Der Untersuchungsbereich 13 hat die Form eines Zylinders, 
der im folgcndcn als Objcktzylindcr 13 bczcichnct wird. 
Nach dem Durchsetzen des Objektzylinders 13 trifft das 60 
Rontgenstrahlenbundel 4 auf eine an der Gantry 1 befestigte 
zweidimensionale Detektoreinheit 16, die eine Anzahl von 
Detektorzeilen mit jeweils einer Vielzahl von Detektorele- 
menten umfaBt. Die Detektorzeilen befinden sich in zur Ro- 
tationsachse senkrechten Ebenen auf einem Kreisbogen um 65 
die Strahlenquelle S; sie konnen aber auch einen Kreisbogen 
um die Rotations achse beschreiben oder geradlinig sein. Je- 
des Detektorelement liefert in jeder Strahlenquellenposition 



einen MeBwert fiir einen Strahl aus dem Strahlenbundel 4. 

Der mit a max bezel chnete Offnungswinkel des Strahlen- 
biindels 4 (als Offnungswinkel ist der Winkel definiert, den 
ein Strahl den einer zur Rotationsachse 14 senkrechten 
Ebene am Rande des Strahlenbiindels 4 liegt, mit einer 
durch die Strahlenquelle S und die Rotationsachse 14 defi- 
nicrtcn Zcntralstrahlcbcnc cinschlicBt) bcstimmt dabci den 
Durchmesser des Objektzylinders 13, innerhalb dessen sich 
das zu untersuchende Objekt bei der Akquisition der MeB- 
werte befindet. Der Untersuchungsbereich 13 - bzw. das 
Objekt oder der Patientenlagerungstisch - kann mittels eines 
Motors 5 parallel zur Rotationsachse 14 bzw. zur z-Achse 
verschoben werden. Die Geschwindigkeit dieses Vorschubs 
in z-Richtung ist vorzugsweise konstant und einstellbar. 
Wenn die Motoren 5 und 2 gleichzeitig laufen, ergibt sich 
eine helixformige Abtastbewegung der Strahlenquelle S und 
der Detektoreinheit. Wenn hingegen der Motor 5 fiir den 
Vorschub in z-Richtung stillsteht und der Motor 2 die Gan- 
try einzeln rotieren laBt, ergibt sich eine kreisformige Ab- 
tastbewegung der Strahlenquelle S und der Detektoreinheit 
relativ zum Untersuchungsbereich 13. Im folgenden soil nur 
diese kreisformige Abtastbewegung betrachtet werden. 

Die von der Detektoreinheit akquirierten MeBdaten wer- 
den einem Bildverarbeitungsrechner 10 zugefiihrt, der dar- 
aus die Absorptions verteilung in einem Teil des Untersu- 
chungsbereichs 13 rekonstruiert und z. B. auf einem Moni- 
tor 11 wiedergibt. Die beiden Motoren 2 und 5, der Bildver- 
arbeitungsrechner 10, die Strahlenquelle S und der Transfer 
der MeBdaten von der Detektoreinheit 16 zum Bildverarbei- 
tungsrechner 10 werden von einer geeigneten Kontrollein- 
heit gesteuert. 

Fig. 2 zeigt den Ablauf eines MeB- und Rekonstruktions- 
verfahrens, das mit dem Computertomographen nach Fig. 1 
durchgefiihrt werden kann. 

Nach der Initialisierung im Block 101 rotiert die Gantry 
mit einer konstanten Winkelgeschwindigkeit. Im Schritt 102 
wird die Strahlung der Strahlenquelle S eingeschaltet, und 
die dabei von den Detektorelementen der Detektoreinheit 16 
akquirierten MeBwerte werden in einem Speicher des Bild- 
verarbeitungsrechners 10 gespeichert. 

Fig. 3 zeigt die kreisformige Bahn 17, auf der sich die 
Strahlenquelle S und die Detektoreinheit 16 um die Rotati- 
onsachse 14 bewegen. Fiir eine bestimmte Strahlenquellen- 
position So ist das Strahlenbundel 4 dargestellt. Man kann 
sich dieses kegelformige Strahlenbundel 4 aus einer Viel- 
zahl von ebenen Strahlenfachera zusammensetzen, die sich 
- wie die Strahlenfacher 401. . .403 in zur Rotationsachse 
parallelen Ebenen befinden. Obwohl man die Strahlen des 
kegelformigen Strahlenbiindels 4 auch auf andere Weise zu- 
sammenfassen kann, soli der Begriff "Strahlenfacher" irn 
folgenden nur fur solche Strahlen verwendet werden, die in 
einer gemeinsamen, zur Rotationsachse 14 parallelen Ebene 
liegen. Diese Strahlenfacher gehen von einer Strahlenquel- 
lenposition aus und werden von je einer zur Rotationsrich- 
tung parallelen Spalte von Detektorelementen auf der De- 
tektoreinheit 16 erfaBt. Fig. 3 deutet an, daB auch in anderen 
Positionen der Strahlenquelle (z. B. S_i, Si oder S2) das 
emittierte kegelformige Strahlenbundel gemessen wird. 

Im Schritt 103 wird der Durchmesser des darzustcllcndcn 
Bereichs FOV (field of view) festgelegt. Der Durchmesser 
kann dem durch <x max definierten Durchmesser des Untersu- 
chungsbereichs 13 entsprechen; er kann jedoch auch kleiner 
sein. 

In den Schritten 104-107 erfolgt ein Rebinning. Die Da- 
ten werden dabei so umsortiert und uminterpoliert, als wa- 
ren sie mit einer anderen Strahlenquelle (einer kreisbogen- 
formigen Strahlenquelle, die zueinander parallele Strahlen- 
facher emittiert) und mit einem anderen Detektor (einem 
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ebenen, reehteckigen und die Rotations achse enthaltenden 
'Virtuellen'' Detektor) gemessen worden. 

Im Schritt 104 werden dazu zunachst die Strahlenfacher 
von verschiedenen Strahlenquellenpositionen zu jeweils ei- 
ner Gruppe zusammengefaBt, die in zueinander parallelen 
Ebenen liegen. Die zu einer Gruppe gehorenden Strahlenfa- 
cher crfullcn dahcr die Bcdingung, daB die Summc dcr Win- 
kel a und p fur alle Strahlenfacher dieser Gruppe den glei- 
chen Wert hat. Dabei ist a der Winkel, den die Ebene des 
Strahlenfachers mit einer durch die Strahlenquellenposition 
und die Rotations achse definierten Zentralstrahlebene ein- 
schlieBt und der durch die Lage der Spalte von Detektorele- 
menten gegeben ist, die den betreffenden Strahlenfacher er- 
faJ3t hat. (3 ist ein die Strahlenquellenposition (z. B. So) auf 
dem Kreis 17 kennzeichnender Winkel. Wenn die Strahlen- 
facher fur eine Strahlenquellenposition diese Bedingung 
nicht exakt erfullen, muB fiir diese Strahlenquellenposition 
ein entsprechender Strahlenfacher durch Interpolation aus 
den Strahlen benachbarter Strahlenfacher ermittelt werden. 

Fig, 4 zeigt eine auf diese Weise gebildete Gruppe von 
Strahlenfachern. Von jeder der einander benachbarten Strah- 
lenquellenpositionen S 2- • -So- • -$2 gehort jeweils ein 
Strahlenfacher zu einer Gruppe. Jede Gruppe kann durch 
den Winkel 3 der Strahlenquellenposition (z. B. S 0 ) gekenn- 
zeichnet werden, deren zu der Gruppe gehorender Strahlen- 
facher die Rotationsachse 14 durchsetzt. (In der Regel ist 
dies die mittlere Strahlenquellenposition - im Beispiel nach 
Fig. 4 also die Strahlenquellenposition S 0 ). Es konnen sich 
dann so viele Gruppen von Strahlenfachern ergeben, wie es 
Strahlenquellenpositionen gibt, aber auch mehr oder weni- 
ger. 

Die auf diese Weise bestimmten Strahlenfacher - darun- 
ter die in Fig. 4 dargestellten Strahlenfacher 411. . .415 - de- 
finieren ein Strahlenbundel 410, das eine zeltartige Form hat. 
und sich aus Strahlenfachern zusammensetzt, die in zuein- 
ander und zur Rotationsachse parallelen Ebenen liegen. In 
Fig. 4 ist auBerdem der Schnittbereich 420 dargestellt, der 
sich ergibt, wenn das Strahlenbundel 410 von einer die Ro- 
tationsachse 14 enthaltenden und zu den Ebenen der Strah- 
lenfacher 411. . .415 senkrechten Ebene geschnitten wird. 
Wie insbesondere Fig. 5 zeigt, die diesen Schnittbereich 
darstellt, sind die oberen und unteren Rander gewolbt. Diese 
Wolbung ist darauf zuruckzufuhren, daB die Strahlenquel- 
lenpositionen in der Mitte (z. B. So) weiter von der Schnitt- 
ebene entfernt sind als die am Rande S2 oder S_2 und daB die 
Facher alle den gleichen Offnungswinkel haben, weil die 
Detektorzeilen einen Kreisbogen um die Strahlenquelle S 
beschreiben. Bei einer anderen Geometrie der Detektorzei- 
len wiirde sich die Form der Schnittebene 420 andern. Bei 
einer ebenen Detektoreinheit (z. B. mit gradlinigen Detek- 
torzeilen) ware die Wolbung noch ausgepragter, weil die am 
Rande liegenden Strahlenfacher, wie z. B. 411 und 415, da- 
bei einen geringeren Offnungswinkel hatten. 

Fiir jede Gruppe von Strahlenfachern wird daher in dem 
ebenen Schnittbereich 420 ein rechteckiger virtueller Detek- 
tor definiert (Schritt 105), dessen oberer und unterer Rand 
161 bzw. 162 durch die Abmessungen der auBeren Strahlen- 
facher 411 bzw. 415 in den ebenen Schnittbereich gegeben 
sind. Wenn dcr Durchmcsscr des zu rckonstruicrcndcn FOV 
(field of view) kleiner gewahlt worden ist als der Durchmes- 
ser des Untersuchungsbereichs 13, konnen von jeder Gruppe 
jeweils die auBeren Strahlenfacher (z. B. 411, 415) entfallen. 
Der Abstand der oberen und unteren Rander 161 und 162 
kann in diesem Fall groBer gewahlt werden als bei dem in 
Fig. 5 dargestellten Beispiel. 

In Fig. 5 sind auBerdem - durch runde Punkte markiert - 
die DurchstoBpunkte einiger in den Strahlenfachern 
411. . .415 enthaltener Strahlen durch diesen virtuellen De- 



tektor dargestellt. SchlieBlich sind durch Kreuze die Stiitz- 
punkte eines regelmaBigen kartesischen Gitters angedeutet. 
Die DurchstoBpunkte und die Stiitzpunkte fallen durchweg 
nicht zusammen. Man erkennt einerseits, daB die Strahlenfa- 
5 cher auBen enger beieinanderliegen als innen und daB die 
DurchstoBpunkte eines Strahlenfachers innen weiter von- 
cinandcr entfernt sind als auBen. Aus den McBwcrtcn fiir die 
DurchstoBpunkte miissen daher in den beiden folgenden 
Schritten 106 und 107 die MeBwerte an den aquidistanten 
Stutzstellen innerhalb des virtuellen Detektors 160 ermittelt 
werden. 

Im Schritt 106 erfolgt dabei zunachst eine vertikale Inter- 
polation derart, daB die Stutzstellen in vertikaler Richtung 
fiir alle Strahlenfacher voneinander den gleichen Abstand 
haben, wie die DurchstoBpunkte bzw. Stutzstellen am Rande 
des virtuellen Detektors 160. Im Schritt 107 erfolgt eine In- 
terpolation in horizontaler Richtung, so daB man interpo- 
lierte Werte fiir Stutzstellen erhalt, die innerhalb des virtuel- 
len Detektors 160 in horizontaler Richtung voneinander den 
gleichen Abstand haben. Die Interpolation sschritte 105 und 
106 konnen auch miteinander vertauscht und ggf. sogar mit- 
einander kombiniert werden. 

Das Ergebnis dieser Interpolation in horizontaler und ver- 
tikaler Richtung ist, daB innerhalb des virtuellen Detektors 
160 die Intensitat der Strahlung an den aquidistanten Stiitz- 
punkten eines regelmaBigen kartesischen Gitters vorliegt. 
Die Stutzstellen definieren ihrerseits wiederum neue Strah- 
lenfacher, die in parallelen Ebenen mit jeweils dem gleichen 
Abstand voneinander liegen. Nur diese neuen Strahlenfa- 
cher (die teilweise mit den origin alen Strahlenfachern iden- 
tisch sein konnen) werden fur die weitere Rekonstruktion 
benutzt. 

Ein Teil der Strahlen, deren DurchstoBpunkte auBerhalb 
des virtuellen Fensters liegen, werden in den Schritten 106 
und 107 nicht benutzt. Es ist daher von Vorteil, den Kollima- 
tor 3 von vornherein so zu gestalten, daB das kegelformige 
Strahlenbundel von vornherein nicht enthalt. Dadurch 
wurde die Strahlenbelastung fiir den Patienten verringert. 
Anstelle von geradlinigen senkrecht zur Rotationsachse ver- 
laufenden Kanten muBte die Kollimatoranordnung zu die- 
sem Zweck nach innen gewolbte Kanten aufweisen, so daB 
die Strahlenfacher, die die Rotationsachse schneiden oder 
davon einen geringen Abstand haben, einen geringeren Off- 
nungswinkel hatten, als Strahlenfacher am auBeren Rand 
des Strahlenbundels. 

Nach den Schritten 103 bis 107 ist somit fiir jede Gruppe 
von Strahlenfachern die Strahlenintensitat an den regelma- 
Bigen Stutzstellen des zu dieser Gruppe gehorenden virtuel- 
len Detektors 160 bestimmt. Dies erleichtert die erforderli- 
che HochpaBfilterung fundamental, weil lediglich eine ein- 
dimensionale Filterung der durch das Rebinning auf den vir- 
tuellen Detektor 160 erzeugten Daten notwendig ist. Im 
Schritt 108 wird daher auf diese Daten eine eindimensionale 
Filterung mit einem rampenformig mit der Frequenz anstei- 
genden Ubertragungsfaktor angewandt. Dazu werden je- 
weils lediglich die Werte von in horizontaler Richtung auf- 
einander folgenden Stutzstellen herangezogen. Diese Filte- 
rung wird fiir alle Gruppen von (neuen) Strahlenfachern 
ausgefuhrt. 

Die nach dem Rebinning und der Filterung fiir die durch 
die Stutzstellen im virtuellen Fenster 160 definierten Strah- 
len ermittelten Daten werden anschlieBend zur Rekonstruk- 
tion der Absorptionsverteilung im Untersuchungsbereich 
durch eine Riickprojektion herangezogen. Durch das voran- 
gegangenen Schritte 103-108 ist es dabei aber - im Gegen- 
satz zu dem eingangs erwahnten bekannten Verfahren - 
nicht erforderlich, die Beitrage einzelner Strahlen zu dem 
Ab sorptions wert eines Voxels mit einem Faktor zu gewich- 
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ten werden, der vom Ab stand des betreffenden Voxels von 
der Strahlenquellenposition abhangt. 

Im Prinzip wurde es fiir die Riickprojektion geniigen, 
wenn man fiir jedes Voxel Daten aus Gruppen heranzieht, 
deren Strahlenfacher miteinander einen Winkelbereich von 5 
insgesamt nur 180° definieren. Es gibt jedoch im Untersu- 
chungsbcrcich cine Viclzahl von Voxcln, die in samtlichcn 
Strahlenquellenpositionen von Strahlung getrofFen worden 
sind und es ist empfehlenswert fiir diese Voxel samtliche 
verfiigbaren MeBwerte heranzuziehen, um ein moglichst gu- 10 
tes Signal/Rauschverhaltnis zu erreichen. Deshalb wird im 
folgenden danach differenziert, ob die Voxel aus alien Strah- 
lenquellenpositionen bestrahlt worden sind oder nicht. 

Nach der Vorgabe eines Voxels im Untersuchungsbereich, 
dessen Absorption rekonstruiert werden soli (Schritt 109), 15 
wird im Schritt 110 zunachst gepriift, ob das betreffende Vo- 
xel aus alien Strahlenquellenpositionen bestrahlt worden ist. 

Dazu wird zunachst auf Fig. 6 hingewiesen, die die durch 
einen Punkt angedeutete Strahlenquelle S, die durch eine Li- 
nie symboiisierte Detektoreinheit 16 und das Strahlenbundei 20 
4 in einer ersten Position bezuglich der Rotationsachse 14 
zeigt. AuBerdem zeigt Fig, 6 in einer zweiten, gegeniiber der 
ersten um 180° versetzten Position die Strahlenquelle S', die 
Detektoreinheit 16' und das Strahlenbundei 4'. AuBerdem ist 
in Fig. 6 der Untersuchungsbereich 13 angedeutet. Die 25 
Strahlen an den in Richtung der Rotationsachse 14 gegen- 
einander versetzten Randern des Strahlenbundels 4 sind ge- 
strichelt angedeutet. Mit ausgezogenen Linien sind hinge- 
gen die Strahlen 40 bzw. 40' angedeutet, die nach dem Weg- 
lassen der Strahlen am Rande (Schritt 107) das Strahlenbun- 30 
del begrenzen. 

Es laBt sich zeigen, daB alle Voxel innerhalb einer ebenen 
kreisformigen Scheibe Vo aus einem Winkelbereich von zu- 
mindest 180° bestrahlt worden sind. Die Seitenflachen die- 
ser Scheibe sind durch die Schnittpunkte der Randstrahlen 35 
40 bzw. 40' mit der Rotationsachse 14 definiert. Alle Voxel, 
die innerhalb des durch die Randstrahlen 40 bzw. 40' defi- 
nierten, diskusformigen Teilvolumens Vi liegen, sind dabei 
aus samtlichen Strahlenquellenpositionen bestrahlt worden, 
wahrend die Voxel, die in dem durch die Differenz von V 0 40 
und Vi gegebenen Teilvolumen V 2 liegen, die Strahlen- 
quelle bei der Abtastbewegung aus einem Winkel von min- 
destens 180° aber weniger als 360° "gesehen" haben. 

Ergibt nun die Abfrage 110, daB das im Schritt 109 ausge- 
wahlte Voxel im Teilvolumen Vi liegt, daB es also aus alien 45 
Strahlenquellenpositionen bestrahlt worden ist, dann erfolgt 
im Schritt 111 eine Riickprojektion, wobei von alien Grup- 
pen jeweils derjenige Strahl in dieser Gruppe herangezogen 
wird, der durch das betreffende Voxel verlauft. Passiert kein 
Strahl exakt die Mitte dieses Voxels, dann muB der zugeho- 50 
rige Wert durch Interpolation der MeBwerte von benachbar- 
ten Strahlen ermittelt werden. 

Ist im Schritt 110 die Abfrage negativ, dann erfolgt im 
Schritt 112 eine weitere Abfrage, ob das ausgewahlte Voxel 
im Teilvolumen V 2 liegt (ob es also die Strahlenquelle unter 55 
einem Winkel von wenigstens 180° "gesehen" hat) ist dies 
nicht der Fall, dann verzweigt das Programm zu einer weite- 
ren Abfrage 113, in der festgestellt wird, ob die Absorption 
schon in samtlichcn Voxcln des FOV rekonstruiert worden 
ist. Ist dies nicht der Fall, dann werden die Schritte 109 ff er- 60 
neut durchlaufen. Bei einer positiven Antwort, d. h. wenn 
sich das im Schritt 109 ausgewahlte Voxel in dem Volumen 
V 2 befindet, werden im Schritt 118 die Gruppen ausgewahlt, 
die zur Rekonstruktion der Absorption in dem Voxel heran- 
gezogen werden. Dazu wird auf Fig. 7 verwiesen, die die 65 
kreisformige Abtastbahn 17 zeigt und auf ihr- durch Punkte 
angedeutet - die verschiedenen Strahlenquellenpositionen. 
Die in Fig. 4 mit S_2. . .Sq. . .S2 bezeichneten Strahlenquel- 
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lenpositionen sind in Fig. 7 durch die zugehorigen Winkel 
P_2. . .po. . .p2 auf dem Kreis 17 gekennzeichnet. Wie in 
Verbindung mit den Schritten 103 bis 108 erlautert, ist jeder 
Strahlenquellenposition eine Gruppe von durch das Rebin- 
ning entstandenen gefilterten MeBdaten zugeordnet. 

AuBerdem ist in Fig. 7 ein Voxel P 2 dargestellt, das - wie 
in Fig. 6 angedeutet - in dem Teilvolumen V 2 liegt. Mit E 2 
ist dabei die zur Zeichenebene der Fig. 7 senkrechte Ebene 
bezeichnet, die durch das Voxel P 2 und die Rotationsachse 
14 definiert wird. SchlieBlich ist mit E eine Ebene bezeich- 
net, die die Ebene E 2 in der Rotationsachse schneidet. 

Es laBt sich zeigen, daB das Voxel P 2 aus alien Strahlen- 
quellenpositionen jenseits der Ebene E bestrahlt worden ist. 
Hingegen hat das Voxel P 2 allenfalls von einigen - nicht 
aber von alien - diesseits der Ebene E (in bezug auf das Vo- 
xel) befindlichen Strahlenquellenpositionen Strahlung emp- 
fangen. 

Infolgedessen werden im Schritt 119 alle durch das Re- 
binning und das Filtern entstandenen Gruppen von Daten 
zur Rekonstruktion herangezogen, die den jenseits der 
Ebene E befindlichen Strahlenquellenpositionen (z. B. 
P2. . .p 2) zugeordnet sind und die miteinander einen Be- 
strahlungswinkelbereich von genau 180° bilden. Dieser Be- 
strahlungswinkelbereich ist notwendig und hinreichend, um 
die Absorption in dem betreffenden Voxel zu rekonstruie- 
ren, auch wenn die mit diesen Daten durchgefuhrte Riick- 
projektion ein schlechteres Signal/Rauschverhaltnis ergibt 
als die im Schritt 111 durchgefuhrte Riickprojektion. 

Sobald im Schritt 113 festgestellt ist, daB alle Voxel in 
dem ausgewahlten FOV rekonstruiert worden sind, ist das 
Verfahren damit beendet (Schritt 114). Die Absorptions- 
werte in dem FOV konnen dann - ggf. scheibenweise - auf 
einem geeigneten Display wiedergegeben werden. 

Patentanspriiche 

1. Computertomographie- Verfahren mit den Schritten 

- Erzeugen eines kegelformigen, einen Untersu- 
chungsbereich bzw. ein darin befindliches Objekt 
durchsetzenden Strahlenbundels mit einer Strah- 
lenquelle, 

- Erzeugung einer kreisformigen, eine Rotation 
um eine Rotationsachse umfassenden Relativbe- 
wegung zwischen der Strahlenquelle einerseits 
und dem Untersuchungsbereich bzw. dem Objekt 
andrerseits, 

- Akquisition von MeBdaten, die von der Intensi- 
ty in dem Strahlenbiindel jenseits des Untersu- 
chungsbereiches abhangen, mit einer Detektorein- 
heit wahrend der Relativbewegung, 

Rekonstruktion der raumlichen Verteilung der 
Absorption innerhalb des Untersuchungsberei- 
ches aus den akquirierten MeBdaten 
gekennzeichnet durch folgende Schritte: 

Rebinning der MeBdaten zu einer Anzahl von 
Gruppen, wobei jede Gruppe mehrere zueinander 
und zur Rotationsachse parallele Ebenen umfaBt, 
in denen sich je ein Strahlenfacher befindet, 
Eindimcnsionalc Filtcrung der durch das Rebin- 
ning erzeugten Daten einer jeden Gruppe in Rich- 
tung senkrecht zur Richtung der Ebenen, 
Rekonstruktion der raumlichen Verteilung der 
Absorption durch Riickprojektion der gefilterten 
Daten von mehreren Gruppen. 

2. Computertomographie- Verfahren nach Anspruch 1, 
dadurch gekennzeichnet daB das Rebinning auf je ei- 
nen zu den Ebenen einer jeden Gruppe senkrechten, 
ebenen und rechteckigen virtuellen Detektor erfolgt, 
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der die Rotations achse enthalt. 

3. Computertomographie-Verfahren nach Anspruch 1, 
gekennzeichnet durch die Rekonstruktion der Absorp- 
tion der Voxel mit einem Bestrahlungswinkelbereich 
von mindestens 180° unter Heranziehung der gefilter- 5 
ten MeBdaten von ausschlieBlich denjenigen Strahlen- 
qucllcnpositioncn, die von dem Voxel durch cine die 
Rotations achse enthaltende Ebene getrennt sind, wel- 
che zu einer durch das betreffende Voxel und die Rota- 
tionsachse definierten Ebene senkrecht steht. 10 

4. Computertomographie-Verfahren nach Anspruch 1, 
gekennzeichnet durch folgende Schritte: 

a) Zuordnung von Voxeln mit einem Bestrah- 
lungswinkelbereich von 360° zu einem ersten 
Teilvolumen, 15 

b) Zuordnung von Voxeln mit einem Bestrah- 
lungswinkelbereich von mindestens 180° aber 
weniger als 360° zu einem zweiten Teilvolumen, 

c) Rekonstruktion der Absorption der Voxel in 
dem ersten Teilvolumen unter Heranziehung von 20 
MeBdaten von samtliche Strahlenquellenpositio- 
nen, 

d) Rekonstruktion der Absorption der Voxel in 
dem zweiten Teilvolumen unter Heranziehung der 
gefilterten MeBdaten von gerade so vielen Strah- 25 
lenquellenpositionen,, daB sich ein Bestrahlungs- 
winkelbereich von 180° ergibt. 

5. Computertomograph zur Durchfiihrung des Verfah- 
rens nach Anspruch s 1 mit einer Strahlenquelle und ei- 
ner damit gekoppelten Detektoreinheit sowie mit einer 30 
Antriebsanordnung, um ein im Untersuchungsbereich 
enthaltenes Objekt und die relativ zueinander um eine 
Rotations achse rotieren zu lassen und mit einer Rekon- 
struktionseinheit zur Rekonstruktion der raumlichen 
Verteilung der Absorption innerhalb des Unter su- 35 
chungsbereiches aus den von der Detektoreinheit ak- 
quirierten MeBdaten, 

gekennzeichnet durch 

a) Mittel zum Re binning der MeBdaten zu mehre- 
ren Gruppen, wobei jede Gruppe mehrere zur Ro- 40 
tationsachse parallele Ebenen umfaBt, in denen 
sich je ein Strahlenfacher befindet, 

b) Mittel zur eindimensionalen Filterung der 
durch das Rebinning erzeugten Daten einer jeden 
Gruppe in Richtung senkrecht zu zur Richtung 45 
den Ebenen 

c) Mittel zur Rekonstruktion der raumlichen Ver- 
teilung der Absorption aus den gefilterten Daten 
verse hiedener Gruppen 

6. Computertomograph, insbesondere nach Anspruch 50 
5, dadurch gekennzeichnet, daB die in Richtung der 
Rotations achse gegeneinander versetzten Rander der 
Kollirnatoranordnung so geformt sind, daB das kegel- 
formige Strahlenbundel in seiner Mitte eine geringere 
Offnung hat als an seinen Randern. 55 
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